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Abstract 



We report trigonometric parallaxes for water masers in the G012. 88+0.48 region and in the massive star forming complex W33 
(containing G012.68~0.18, G012.81~0.19, G012.90-0.24, G012.90-0.26), from the Bar and Spiral Structure Legacy (BeSSeL) survey 
using the Very Long Baseline Array. The parallax distances to all these masers are consistent with 2.40^55"]j kpc, which locates the 
W33 complex and G012. 88+0.48 in the Scutum spiral arm. Our results show that W33 is a single star forming complex at about 
two-thirds the kinematic distance of 3.7 kpc. The luminosity and mass of this region, based on the kinematic distance, have therefore 
been overestimated by more than a factor of two. The spectral types in the star cluster in W33 Main have to be changed by 1.5 points 
to later types. 
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1. Introduction 

W33 is a massive and very luminous star forming complex that 
contains typical star forming regions from quiescent infrared 
dark clouds to highly active, infrared bright clouds, associated 
with H II regions. Observations of the dust emission in the W33 
complex from the APEX Teles cope Large A rea Survey of the 
GALaxy (ATLASGAL, Schulle r et al.l l2009l) at 870 fim show 
three large molecular clouds and several smaller dust clumps 
(Fig. [TJ- Each of the large clouds contains water masers at 
G012.68-0.18 (W33B), G012.81-0.19 (W 33Main). G012.90- - 
0.24 ( W33 A), and G01 2.90-0.26 (W33 A) dGenzel & Downesl 
[19771 iJaffeetalLflgsTI) . 

To derive fundamental physical parameters of the molecu- 
lar clouds such as their size, luminosity, mass, and the spec- 
tral types of embedded stars, we have to know their distances. 
One commonly used method is to dermine the kinematic dis- 
tance of a cloud from radial velocity measurements. Comparing 
the radial velocity of the cloud with a rotation curve of the 
Galaxy (e.g. Burton & Gordon, 1978; Reid et al., 2009b) yields 
the Galactocentric radius from which the distance of the cloud 
can be determined. However, in the first and fourth quadrants 
there are two possible distances corresponding to a given ra- 
dius, complicating the distance determination. In addition, local 
velocity deviations due to shocks, outflows or noncir cular mo- 
tions can lead to large errors of the kinematic distance (IXu et al.L 
l2006h . 

Kinematic distance determination for the W33 complex is 
complicated by a peculiar kinematic structure that was discov- 
ered via H ?CO absorption line and radio recombination line ob- 
serva tions dBieging et al.L 119781; iGoss et al.L 119781; iBieging et all 
|1982 |). Two radial velocity components at 36 and 58 km s ' were 
detected over extended paits of the region. The first component 



is concentrated in the northeastern part (W33 A) and the sec- 
ond component peaks in the southwestern part (W33 B). This 
velocity structure can be explained by either one connected star 
forming region with large internal motions at a near kinematic 
distance of 3.7 kpc (corresponding to the radial velocity of 36 
km s '), or a superposi tion of independent star forming regions 
along the line of sight (Has chick & 1101119831: IGoss et ali[l978t 
Gardner etaLlflgT S; Biegin g et al.lll978l) . 

Assuming a near kinematic distance of 3.7 kpc to W33, 
the angular size of the complex of 15' c an be converted into a 
physical size of 16 pc. lStier et al.l d 19841) determined a total in- 
frared luminosity of ~2-10^ Lq for W33 at this distance. Using 
'^C O maps and a typica l value of n('^CO)/n(H2) = 4 • 10"\ 
Gol dsmith & Mad (Il983h estimated the total mass of the com- 
plex to be be tween 0.2 and 2-10^ M (T|. From radio continuum 
observations. lHaschick & Ha (Il983h inferred the existence of an 
OB star cluster in W33 Main with individual spectral types rang- 
ing from 06 to BO, based on a kinematic distance of 4 kpc. 

W33A is a well-studied star forming regi on. The source 
has a high far-infrared luminosity of 10'' Lq (iFaundez et al.L 



20041) at the previous l y assu med kinematic distance of 3.7 kpc. 
Galvan-Madrid et aTl (1201 Ol) detected parsec-scale filaments of 
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cold molecular gas around two dust cores. The interaction of 
these filaments might have triggered the star formation activity 
in the cores. The brightest core drives a strong outflow and the 
dynamics in this core show evidence for a rotating disk, perpen- 
dicular to the outflow. These r esults are supported by CO obser- 
vations of iDavies et al.l (12010 ') which also suggest the presence 
of a rotationally-flattened cool molecular envelope. Bry emis- 
sion in W33 A seems to trace a fast bipolar wind which has th e 
same orientation as the large-scale outflow (Davies et al., 120101) . 
The colder and more massive second core in W33 A seems to be 
in an earlier evoluti onary stage, driving a more modest outflow 
(iGalvan-Madrid et a l., 2010). 

G012. 88+0.48 is a high-mass p rotostellar object 
(ISridharan et al.L l2002t iBeuther et idl l2002l) which is prominent 
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Figure 1: 870 //m dust emission of the W33 complex from the ATLASGAL survey with contour levels of 1, 4, 7, 10, 13, 16, 20, and 
30 Jy beam ^ The positions of the water masers are indicated with black crosses. The arrows show the internal motion vectors of 
the two clouds W33 A and W33 B in the reference frame of W33 Main. The length of the an^ow in the upper left corner corresponds 
to a motion vector of 1 mas yr '. The circle in the lower left comer shows the beam of the ATLASGAL observations. G012.88+0.48 
is located at an angular distance of about 0.7 degrees to the north of W33, outside of this figure. 



for its strong maser emission from the OH, H2O, and CH3OH 
mol ecules. G012.88+0.48 has a radial velocity of 33.8 km 
s~' dSridharan et al.L l2002h . locating the star forming region 
at a distance of 3.6 kpc, assuming it is at its near-kinematic 
distance. This soure is projected about 0.7 degrees from the 
main star forming re gions of W33. Sridharanet al. (2002) and 
iBeuther et al.l (l2002l) determined the luminosity and mass of 
this star forming region to be L ~ 3-10^ L0 and M ~ 310-' Mq 
at this distance. IXu et alJ (l201lh determined the trigonometric 
parallax for the 12.2 GHz methanol masers in this source. They 
obtained a value of 0.428+0.022 mas which corresponds to a 
distance of 2.34 kpc, revising the kinematic distance by a factor 
of 0.65 (= 2.34kpc/3.6 kpc). If G012.88+0.48 is associated with 
the W33 complex, then the kinematic distances of the complex 
are too large. Alternatively, it is possible that G012.88+0.48 is 
not associated with W33 and is substantially closer to the Sun. 

Other recent trigonometric parallax measurements have 
shown that kinematic distance measurements can be wrong by 
more than a factor of two (e.g. Reid et al., 2009b; Sato et al., 
r2010al) . Trigonometric parallax observations of maser sources, 
however, have proven to be a reliable method to accurately mea- 
sure distances to star forming regions whose accuracy is limited 
only by observational uncertainties. 

The aim of this project is to obtain distances to the water 
masers in the W33 complex from trigonometric parallax ob- 
servations and, therefore, address the question of whether the 



molecular clouds in W33 belong spatially to one connected star 
forming complex or are at different distances and projected near 
each other on the sky. Also, solid distances will provide the ba- 
sis for more accurate values for the luminosity and mass of the 
complex. 

The observations presented in this paper are part of 
the Bar and S pi ral St ructure Legacy (BeSSeL) Survejo 
jBrunthaler et al.L |2011|) . a National Radio Astronomy 
Observatory'^ Very Long Baseline Array (VLBA) key 
project which will determine the distances to hundreds of star 
forming regions. We report results of the parallax and proper 
motion measurements toward the water masers in the sources 
G012. 88+0.48 and the massive star forming complex W33. 



2. Observation and Data Reduction 

The massive star forming region G012.88+0.48 and the 
W33 complex (G012.68-0.18, G012.81-0.19, G012.90-0.24, 
G012.90-0.26) were observed in the 616 - 523 transition of the 
H2O molecule (rest frequency 22.23508 GHz), using the VLBA. 
The general observation setup and data calibration procedures 
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are described in lReid et alJ (l2009al) . Here, we only describe the 
details specific to the present sources. 

All water masers were observed at nine epochs (VLBA pro- 
gram BR 1451) with a total observing time per epoch of ~ 7 
hours. The observing dates (2010 Sep 19, 2010 Oct 22, 2010 Dec 
19, 201 1 Feb 15, 201 1 Mar 26, 201 1 May 21, 201 1 Jul 30, 201 1 
Sep 20, 2012 Jan 03) were selected to well sample the Right 
Ascension parallax signature in time, since the Declination par- 
allax signature is much smaller. However, due to bad weather 
conditions, the data of the first epoch could not be used for the 
parallax measurements. 

Three diff^erent quasars, J1808-1822, J1809-1520, and 
J 1825- 17 18, from the VLB A caHbrator search for th e BeSSeL 
surve y dimmer et al.L 1201 Ih and the ICRF2 catalog dPev et al.L 
[2009), were observed for background position references. The 
background source J1809-1520 could only be detected in the 
first four epochs, and thus was not used for the parallax determi- 
nation. Four adjacent frequency bands with 4 MHz bandwidth 
each were used in right and left circular polarization. The maser 
signals were centered in the second band at an LSR velocity, 
Vlsr, of 45 km s"' for all masers. The channel spacing of the 
observations is 0.42 km s"'. 

In order to correct for instrumental phase offsets between 
the frequency bands, observations of the calibrator J1800+3848 
from the ICRF2 catalog were used. The strongest maser feature 
in the water maser G012.68-0.18 at a Wlsr of 59.3 km s"' was 
used for the phase-referencing of the data. 

The water masers G012.81-0.19, G012.90-0.24, G012.90- 
0.26, and G012.88+0.48 were much weaker than G012.68-0.18 
and were calibrated using G012.68-0.18 as phase reference. 
This yields positions relative to G012.68-0.18, from which rela- 
tive parallaxes can be estimated. 

For the da ta reduction, we used the software Parseltongue 
dKettenis et al.. 2006), a scripting interface for NRAO's 
Astronomical Image Processing System (AIPS). After calibrat- 
ing the data, the maser emission and the continuum sources were 
imaged with the AIPS task IMAGR with a circular Clean beam 
of 2 mas for all epochs. The positions of the masers and the 
background sources J1808-1822 and J1825-1718 were then de- 
termined by fitting Gaussia n brightness distributions to the im- 
ages (see lReid et al.Ll2009al) . Absolute positions of the strongest 
maser spot in G 12.68-0. 81 were derived relative to the two back- 
ground sources. For the remaining water masers, relative posi- 
tions to G012.68-0.18 were determined. To obtain absolute par- 
allaxes, we added the positions of G012.68-0.18 relative to both 
background quasars to the position measurements of these water 
masers. 

The positions were modeled with a sinusoidal parallax sig- 
nature and linear proper motions in each coordinate. Since sys- 
tematic errors (typically from unmodeled atmospheric delays) 
are normally much larger than random noise, the formal posi- 
tion eiTors from the Gaussian fits are mostly too small to ac- 
count for the total error of the observations which leads to high 
X^ values in the modeling process. The magnitude of the system- 
atic errors can only be inferred from the quality of the parallax 
fit. Thus, we defined "error-floors" for the Right Ascension and 
Declination maser positions and added those in quadrature to the 
Gaussian fit errors. The "eiTor-floor" values were adjusted until 
the x^ value per degree of freedom was close to unity for both 
coordinates (see lReid et aI.Ll20()9a) . 

In order to determine the internal gas kinematics in each 
source, we measured the proper motions for all maser spots that 
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were detected in at least three epochs without solving for the 
parallax parameter. 



3. Results 

Table[T]gives the coordinates (Cols. 2, 3), line-of-sight velocities 
(Col. 4), and peak flux densities from epoch 201 1 May 21 (Col. 
5) of the background quasars and the maser spots that were used 
for the parallax fits. Column 6 lists the rms brightness of the 
reference channel images, obtained from epoch 2011 May 21. 
In the last column, we show the associations of the masers in the 
W33 complex. 

The parallax and proper motion results of all masers are sum- 
marized in Table |2] Columns 2, 3, and 4 give the line-of-sight 
velocities of the masing and thermal gas and the parallaxes of 
the maser sources. In the sixth and seventh columns, the proper 
motions in eastward and northward directions are listed. 

One goal of our observations is to estimate the distance and 
proper motion of the central object in each source. In G012.68- 
0. 18, we obtained the proper motion o f the central object by fit- 
ting a model for expanding outflows dSato et al.L 1201 Obh to the 
proper motions of all maser spots. In the other sources, the small 
number of maser spots and their distribution did not allow the fit- 
ting of this model to their proper motions. Instead, the presented 
values are averages of the proper motions of all maser spots in 
each source. 



3.1. W33 

Figure [T] shows the dust emission of the W33 complex from the 
ATLASGAL survey at 870 fim. The three most massive clouds 
W33 A, W33 Main, and W33 B are indicated. The water masers 
are marked with black crosses. 

Water maser spots were detected in an 0."15 x 0."15 area 
in W33 B (G012.68-0. 18), covering a velocity range of 57 - 63 
km s '. Figure |2] shows that the maser emission in the reference 
channel (Vlsr = 59.3 km s"', left panel) as well as the emis- 
sion of the three background quasars J 1808-1 822 (light upper 
panel), J1809-1520 (left lower panel), and J1825-1718 (light 
lower panel) are dominated by one compact component. Since 
the quasar J 1 809-1 520 was detected only in the first four epochs, 
we determined the position of the reference maser spot only rel- 
ative to the background quasars J1808-1822 and J1825-1718 as 
a function of time (see Fig. [3]l. For three epochs (2010 Oct 22, 
2010 Dec 19, 201 1 Sep 20), we used a different maser spot at the 
same velocity as position reference since the original reference 
maser spot, used in the other epochs, was too weak to serve as a 
phase reference. 

The positions of the reference maser spot in G012.68-0.18 
relative to background quasars were fitted with the sinusoidal 
parallax signature and a linear proper motion. We obtained par- 
allax estimates of tt = 0.462 + 0.034 mas and n = 0.387 + 0.037 
mas, relative to the background quasars J 1808-1 822 and J 1825- 
1718, respectively. Since both results are consistent within their 
joint uncertainty, we combined the two data sets and repeated 
the fitting, yielding a value of tt = 0.416 + 0.028 mas which 
corresponds to a distance of 2.40^q[^ kpc (see Fig.|3]l. 

We fitted the proper motions of all maser spots in G012.68- 
0. 18, detected in at least three epochs, with a model for expand- 
ing outflows with the position and proper motion of the cen- 
tral driving s ource as free para meters (for more details about 
the model see lSato et al.Ll2010bl) . The proper motion of the cen- 
tral star is fi_, = -1.00 + 0.30 mas yr"' and fiy = -2.85 + 0.29 



K. Immer et al.: Trigonometric Parallaxes of G012.88+0.48 and W33 
Table 1 : Coordinates, line-of-sight velocities, and peak flux densities of the background quasars and the strongest maser spots. 



Source 


R.A. (J2000) 


Dec. (J2000) 


^LSR 


FpraJr 


rms 


Cloud Association 




[hh mm ss.ssss] 


[ dd " "". ] 


[kms-'] 


[Jybeam-'] 


[Jy beam"'] 




J 1808- 1822 


18 08 55.5150 


-18 22 53.383 




0.012 


0.001 




J 1809- 1520 


18 09 10.2094 


-15 20 09.699 




0.022 


0.001 




J1825-1718 


18 25 36.5323 


-17 18 49.848 




0.096 


0.001 




GO 12.68-0. 18 


18 13 54.7457 


-18 0146.588 


59.3 


132.5 


0.001 


W33B 


G012.81-0.19 


18 14 13.8283 


-17 55 21.035 


34.1/- 1.4 


1.0/0.2 


0.018/0.007 


W33Main 


GO 12.90-0.24 


18 14 34.4366 


-17 51 51.891 


34.9 


18.3 


0.081 


W33A 


GO 12.90-0.26 


18 14 39.5714 


-17 52 00.382 


37.0 


3.0 


0.016 


W33A 


G012.88+0.48 


18 1151.4939 


-17 31 28.911 


29.4 


23.3 


0.133 





Notes. Columns 2, 3, and 4 list the coordinates and line-of-sight velocities of the background quasars and the maser spots that were used for 
parallax fitting. The peak flux densities of the background quasars and the strongest maser spots and the rms of the reference channel images, 
obtained from epoch 2011 May 21, are listed in Cols. 5 and 6. The last column shows with which cloud the masers in the W33 complex are 
associated. 

Table 2: Parallax and proper motion results. 



Maser 
Source 



V 



LS Rjnaser 

km s"' 



* LSRjhermal 

km s"' 



Parallax n 
[mas] 



Distance d 
[kpc] 



/'a 

[mas yr"'] 



[mas yr"'] 



G012.68-0.18(W33B) 


59.3 


-55 


0.416+0.028 


2 40+"" 


-1.00±0.30 


-2.85±0.29 


G012.81-0.19(W33Main) 


-1.4 


-36 


0.343±0.037 


9 Q9+035 
^■^^-0.28 


-0.24±0.17 


+0.54±0.12 


G012.81-0.19(W33Main) 


34.1 


-36 


0.343±0.037 


2 92+"-^^ 


-0.60±0.11 


-0.99±0.13 


G012.90-0.24(W33A) 


34.9 


-36 


0.408±0.025 


9 As+o.ie 


+0.19±0.08 


-2.52±0.32 


G012.90-0.26(W33A) 


37.0 


-37 


0.396±0.032 


2 53+"-^ 


-0.36±0.08 


-2.22±0.13 


G012.88+0.48 


29.4 


-34 


0.340+0.036 


9 04+0.35 


+0.12+0.13 


-2.66+0.23 



Notes. '"' ISridharan et al.l(l2002h : |Purcell et aD(l2012h : IWienen et al.l ( l2012h . Columns 6 and 7 give the proper motions in eastward and northward 
directions. The proper motions of the first maser are obtained from fitting a model of expanding outflows to all maser spots. The proper motions 
of the remaining sources are averages of the proper motions of all maser spots in each source. 



mas yr"'. Figure |4] shows the positions of the maser spots with 
their proper motions in the reference frame of the central source. 
Observations with the Submillimeter Array (SMA) at 230 GHz 
(Immer et al, in prep.) show that the reference maser spot is lo- 
cated in a dust core, offset by ~ 0. 1 " from the submillimeter 
continuum peak. The radial velocity of the reference maser spot 
is similar to the systemic velocit y of the thermal gas in this cloud 
(55 km s-'. IWienenet"ani2012h . 

The remaining W33 masers and G012. 88-1-0. 48 were phase- 
referenced to G12.68-0.18 and the positions of their strongest 
maser spots were determined relative to this maser, yielding rel- 
ative parallaxes. To obtain absolute parallax values, the positions 
of G012.68-0.18 relative to the two background quasars were 
added to the position measurements of the other water masers . 

In G012.81-0.19, we detected two spatially and kinemati- 
cally distinct maser groups, covering velocity ranges of 33 to 35 
km s~' and -5 to km s"', separated by ~21". The positions 
of two maser spots at -1.4 and 34.1 km s"' were determined 
relative to G012.68-0. 18. Combining the two position measure- 
ments, we obtained an absolute parallax of ;7r - 0.343+0.037 mas 
(see Fig. |5j left panel). The positions of the two maser groups 
are shown in Fig. |6] The images on the right-hand side show 
both maser groups on a smaller scale with their proper motions. 
Averaging the proper motions of all maser spots in each group, 
yields a proper motion value for each velocity group (uncor- 
rected for expansion): yU „a^^ „-i = -0.24 ± 0.17 mas yr"' and 

^,34kms-i = 0.54 + 0.12 mas yr"' as well as ^i^_^^^^-^ = 

-0.60 + 0.11 mas yr-' and ju _ji,njs-i = -0.99 + 0.13 mas yr"'. 

The maser spots in the velocity group around -1 km s ' are or- 
dered along a ridge with redshifted velocities at the north-east 



end and blueshifted velocities at the south-west end. The maser 
spot group at 34 km s"' is located in a dust core, offset by ~ 7 " 
from the main continuum emission peak at 230 GHz (Immer et 
al., in prep.). The radial velocity of the reference maser spot is 
similar to the systemic velocity of the thermal gas of 36 km s"' , 
which was estimated from ammonia observations of the H2O 
southern Galactic Plane Survey (HOPS lPurcell et al.Ll20T2l) . The 
maser spot group at -1 km s~' is associated with emission from a 
dust core at 345 GHz (Immer et al., in prep.). However, the radial 
velocities of these maser spots are very different from the sys- 
temic velocity in this source and they might trace the blueshifted 
shock front of an outflow. 

The maser emission from G012.90-0.24 spans a velocity 
range of 33 to 37 km s"' . For the parallax fitting, we used the po- 
sitions of three maser spots at velocities of 34.5, 34.9, and 35.3 
km s"' and obtained an absolute parallax of tt = 0.408 + 0.025 
mas (see Fig. |5] middle panel). The proper motions of the de- 
tected maser spots are very similar (Fig. |7). The proper motion 
of the source (uncorrected for expansion), obtained from aver- 
aging all maser spot proper motions, is ju, = 0.19 + 0.08 mas 
yr"' and yu,, = -2.52 + 0.32 mas yr"'. The ATLASGAL map 
(Fig. [T]) shows that the reference maser spots are offset by ~ 8 
" from a submillimeter continuum peak. The radial velocities of 
the reference maser spots are close to the systemic velocity of 
G012.90-0.24 (36 km s"', estimated from the HOPS survey). 

Only three maser spots were detected in G012.90-0.26, cov- 
ering a velocity range of 36 to 38 km s"' . We used the strongest 
maser spot at 37.0 km s"' for the parallax fitting, yielding an 
absolute parallax of ;7r = 0.396 + 0.032 mas (see Fig. |5] right 
panel). Averaging the proper motions of the three maser spots, 
we obtain a proper motion of the source of //f - -0.36 + 0.08 
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Figure 2: Left upper panel: Maser emission in reference channel 
(59.3 km s '). Right upper, left and right lower panel: Compact 
emission of the three background quasars J 1808-1 822, J 1809- 
1520, and J1825-1718, which were phase-referenced to the ref- 
erence maser channel. The images are obtained from epoch 201 1 
May 21. The contour levels are 10 to 100 % in steps of 10 % of 
the maximum in each panel (G012.68-0.18: 132.5 Jy beam"', 
J1808-1822: 12.3 mJy beam"', J1809-1520: 22.1 mJy beam"', 
J1825-1718: 96.1 mJy beam '). The restoring beam is shown in 
the lower left comer 
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Figure 3: Parallax and proper motion fits of G012.68-0.18, 
yielding a parallax of (0.416+0.028) mas. The positions of 
the reference maser spot at 59.3 (triangles) km s ' were mea- 
sured relative to the background quasars J 1808-1 822 (blue) and 
J1825-1718 (red). Left Panel: Position on the sky with a label 
for the first and last epoch. Middle Panel: East (continuous line) 
and North (dashed line) position offsets with parallax and proper 
motion fits versus time. Right Panel: Right ascension (contin- 
uous line) and Declination (dashed line) parallax fits with the 
best-fit proper motions removed, showing only the parallax sig- 
nature. 
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Figure4: Positions of all maser spots in G012.68-0.18 that have 
been detected in at least three epochs (positions from epoch 201 1 
May 21). The size of the spots scales logarithmically with the 
flux density of the spots. The asterisk marks the position of the 
central star which has an error of ~0.1" in both coordinates. The 
arrows show the proper motions of the maser spots in the refer- 
ence frame of the central star 
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Figure 5: Right ascension parallax fits of G012.81-0.19 (left), 
G012.90-0.24 (middle), and G012.90-0.26 (right), yielding ab- 
solute paraUaxes of 0.343 + 0.037 mas, 0.408 + 0.025 mas, and 
0.396 + 0.032 mas, respectively. The positions of the maser spots 
were measured relative to the water maser G012.68-0.18 and 
consecutively to the two background quasars J 1808-1 822 (blue) 
and J1825-1718 (red). Left panel: Positions of the maser spots 
at 34.1 (triangles) and -1.4 (squares) km s"'. Middle panel: 
Positions of the maser spots at 34.5 (triangles), 34.9 (squares), 
and 35.3 (pentagons) km s '. Right panel: Positions of the maser 
spot at 37.0 (triangles) km s"' . 



mas yr"' and yu,. = -2.22 + 0.13 mas yr"'. The position of the 
reference maser spot is consistent with the submillimeter con- 
tinuum peak MMl-SE, detected by Galvan-M adrid et al. (201^) 
at 230 G Hz. This object is assoc i ated w ith a faint radio source 
at 7 mm (Ivan der Tak & MentenL l2005l) . The radial velocity of 
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G012 81-0 19: (0,0) = 18:14'13. 8^83 -17:55:21.035 (,12000) 



G012,90-0,26; (0,0) = 18:14:39.5714 
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Figure 6: Left panel: Positions of all maser spots in G012.81- 
0.19 that have been detected in at least three epochs (positions 
from epoch 2011 May 21). Inside panels: Zoom-in on the two 
maser spot groups at 34 km s"' (lower panel) and - 1 kms"' (up- 
per panel). The maser spots in the upper panel are ordered along 
a ridge. The size of the spots scales logarithmically with the flux 
density of the spots. The arrows show the absolute proper mo- 
tions of the maser spots. 
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Figure 8: Positions of all maser spots in G012.90-0.26 that have 
been detected in at least three epochs (positions from epoch 201 1 
May 21). The size of the spots scales logarithmically with the 
flux density of the spots. The arrows show the absolute proper 
motions of the maser spots. 



G012.90-0.24; (0,0) = 18:14:34,^ 

\ n 

1 mas/yr 



-17:51:51.891 (J2000) *LSB _ 

km s 



3.2. G012.88+0.48 



;.17 



-o.ooa 

East Offset (arcsec) 



Figure 7: Positions of all maser spots in G012.90-0.24 that have 
been detected in at least three epochs (positions from epoch 201 1 
May 21). The size of the spots scales logarithmically with the 
flux density of the spots. The arrows show the absolute proper 
motions of the maser spots. 



the reference maser spot is very close to th e systemic velocity of 
the th ermal gas in W33A (38.5 km s~'. iGalvan-Madrid et al.L 
l2010h . A comparison of the proper motion of G012.90-0.26 
with CO observations of t he large-scale outflow in W33 A 
(iGalvan-Madrid et al.L l201Cl ) show that their directions do not 
agree. Thus, we will assume that the proper motion of this wa- 
ter maser is not strongly influenced by the outflow but probably 
reflects the motion of the central object. 
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Figure 9: Parallax and proper motion fits of G012.88-t-0.48, 
yielding an absolute parallax of 0.340 ± 0.036 mas. The posi- 
tions of the maser spot at 29.4 (triangles) km s"' were measured 
relative to the water maser G012.68-0.18 and consecutively to 
the two background quasars J1808-1822 (blue) and J1825-1718 
(red). Left Panel: Position on the sky with a label for the first and 
last epoch. Middle Panel: East position offsets with parallax and 
proper motion fits versus time. Right Panel: Right ascension par- 
allax fit with the best-fit proper motions removed, showing only 
the parallax signature. 



Maser emission in the spectrum of G012.88h-0.48 was de- 
tected only in the velocity range 27 to 31 km s"'. Fitting the po- 
sitions of the strongest maser spot at 29.4 km s"', we obtained an 
absolute parallax of :7r = 0.340 ±0.036 mas (see Fig.|9]l, which is 
consistent within 2cr (joint uncertainty) with the parallax result 
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Figure 10: Positions of all maser spots in G012. 88+0.48 that 
have been detected in at least three epochs (positions from epoch 
2011 May 21). The size of the spots scales logarithmically with 
the flux density of the spots. The arrows show the absolute 
proper motions of the maser spots. The smaller picture is a zoom 
on the western group of maser spots. These spots might be trac- 
ing a bow shock in this source. 



of IXu et al] ( l2011h for 12.2 GHz methanol masers in this star 
forming region. 

Figure [To] shows two groups of maser spots, separated by ~ 
0.25". The spatial distribution of the spots in the western group 
shows an arc-like morphology. The combination with the proper 
motions of these spots seems to indicate that the water masers 
trace a shock front in this source. Methanol maser observations 
in G012.88+0.48 at 6.7 GHz with the Japanese VLBI Network 
and the East- Asian VLBI Network present a similar picture with 
a group of maser spots arranged in an arc-like structure with line- 
of-sight velocities from 30 to 40 km s"' (Fuiisawa et al., 2012, 
Fujisawa et al., in prep.). This maser group is offset from our 
western group by -0.3" and -1.2" in east-west and north-south 
direction, respectively. 

We averaged the proper motions of all water maser spots 
to obtain one proper motion for G012. 88-1-0.48, yielding /i^ - 
0.12 + 0.13 mas yr"' and/v,. = -2.66 + 0.23 mas yr"'. 

4. Discussion 

The absolute parallaxes of G012.81-0.19, G012.88+0.48, 
G012.90-0.24, and G012.90-0.26 are consistent with a paral- 
lax of 0.416 mas within 2cr, i.e. the distances to these masers are 
in accordance with the distance to G012.68-0.18 (2.40 kpc). 

Comparing the coordinates of the W33 complex wit h a CO 
latitude-velocity map from the survey of iDame et aP (12001 1) 
(Fig. fTTT i. we find that the complex is associated with a bright 
cloud at the same latitude with a mean velocity of 34.6 km s"', 
which falls along the longitude-velocity locus of the Scutum spi- 
ral arm. 

Since all water masers in the W33 complex are at a sim- 
ilar distance, we conclude that the molecular clouds W33 A, 
W33 Main, and W33 B are connected and belong to the same 
star forming complex. Our parallax distance to this star forming 
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Figure 1 1 : Velocities of the water maser groups in the CO 
latitude-velocity plot (white dots). The CO emission was inte- 
grated over a longitude range from 12.625° to 13.0°. The units 
of the color bar are degrees Kelvin. The water masers in W33 
and in G012. 88-1-0. 48 are associated with clouds at 34.6 km s"' 
and 30.9 km s"', respectively. The water masers at -1 km s"' 
(W33 Main) and 58 km s"' (W33 B) seem to be oudiers. 



cluster is a factor of 0.65 (= 2.4 kpc/3.7 kpc) smaller than the 
kinematic distance. Thus, the luminosity and mass of this com- 
plex have been overestimated in the past by a factor of more than 
two and should be revised to L ~ 8-10^ L© and M ~ (0.8-8)- 10^ 
Mq. 

The spectral types of the stars in the W33 Main cluster were 
determined from the estimated number of Lyman continuum 
photons needed to ionize the material in the surrounding H II re- 
gion. Since the number of photons depends on distance squared, 
the values have been overestimated by a factor of more than two, 
changing the spectral types by 1 .5 points to later types, i.e., a star 
previous estimated to be spectral type 06 should have a spectral 
type 07.5. Thus, the star cluster in W33 Main has spectral types 
ranging from 07.5 to B1.5. 

Our observations show maser detections in the W33 complex 
in three different velocity ranges: 



63kms-i inG012.68-0.18, 

-1 in G012.81-0.19, G012.90-0.24, G012.90- 



57 

33 - 38 km s 

0.26, 

-5-Okms-' 



inG012.81-0.19. 



These findin gs are confirmed by water a nd hydroxyl mase r ob- 
servations bylGenzel & DownesI (Il977h . iLada et all (Il98lh . and 



lArgon et aT (2000). Ammonia observations with the Efifelsberg 
100 meter telescope (Wienen et al., 2012) and from the HOPS 
survey (iPurcell et al.Ll2012h give line-of-sight velocities of ~55, 
-36, -36, and -37 km s"' for G012.68-0.18, G012.81-0.19, 
G012. 90-0.24, and G012.90-0.26, respectively, in accordance 
with our maser observations. On the other hand, lUrquhart et al.l 
(2008) and Chen et al. (2010) detected CO emission peaks at 
35, 60, and 52 km s' in G012.68-0.18, G012.81-0.19, and 
G012.90-0.26. These results show that both velocity compo- 
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nents at -36 and ~58 km s ' are spread over the entire complex 
but the emission at these velocities peaks in different regions. 
The ammonia observations of Wienen et al. (2012) show the ve- 
locity component at ~58 km s"' also at positions in the vicinity 
ofW33 (e.g. G11.94-0.26,G12.74-0.10,G12.90-0.03). This in- 
dicates that this velocity component is not unique to the inner 
W33 complex and probably has a different origin than internal 
motions (e.g. outflows) in the clouds. The velocity component at 
~ -Ikms"' inG012.81-0.19 was not detected in previous ther- 
mal line observations. The arrangement of the water maser spots 
in this velocity group (see Fig. |6]l might suggest that this veloc- 
ity component is associated with the shock front of an outflow in 
the active cloud W33 Main (Immer et al., in prep.) . 

Although water masers are often found at the base of bipolar 
outflows, the small radial velocity difference between the mas- 
ing and the thermal gas in G012.81-0.19, G012.90-0.24, and 
G012.90-0.26 (see Table |2|, the small radial velocity range of 
the masing gas and the proximity of the masers to submillime- 
ter dust peaks suggests that the 3D motions of the water masers 
are similar to the motions of the associated central objects and 
not drastically altered by strong outflows. The only case where a 
strong outflow could affect the proper motion of the water maser 
without changing its radial velocity is if the outflow was aligned 
with the plane of the sky. However, the probability of outflows 
being in the plane of the sky in all three sources is small. Thus, 
we assume that the average proper motion of the water maser 
spots in each source is similar to the motion of its central object. 
Under this assumption, we can estimate the cloud motions inter- 
nal to the W33 complex from the motions of the central objects. 
In G012.68-0.18, we estimated the motion of the central object 
by fitting a model of expanding outflows to the proper motions 
of the maser spots. Since the maser spots in the velocity group 
around -1 km s"' in G012.81-0.19 have a very different ra- 
dial velocity than the thermal gas in this cloud, their proper mo- 
tions probably do not reflect the motion of the central object of 
the W33 Main cloud. They were thus excluded from the follow- 
ing calculations. We averaged the motions of W33 A (G012.90- 
0.24, G012.90-0.26), W33B (G012.68-0.18), and W33Main 
(G012.81-0.19-34kms"'), yielding an estimate of the average 
proper motion of the complex: jij- - -0.44 ± 0.25 mas yr~' and 
//,, - -2.14 + 0.41 mas yr"\ corresponding to -5 km s~' and 
-24 km s ' in the eastward and northward directions, respec- 
tively (at the measured distance of 2.40 kpc). To determine the 
cloud motions of W33 A and W33 B relative to the W33 Main 
cloud, the motion of W33 Main was subtracted from the motion 
of each cloud. Figure [T] shows the cloud motions of W33 A and 
W33 B in the reference frame of W33 Main on the 870 jum dust 
emission map of the W33 complex. 

That W33 A and W33 Main have similar radial velocities 
suggests that most of the motion of W33 A in the reference frame 
of W33 Main is in the plane of the sky, tangentially to W33 Main 
(total speed ^ 17 km s"'). 

The difference in radial velocity between W33 Main and 
W33 B is ~ 22 km s"' . This large velocity component along the 
line of sight coincidentially is equal to the velocity component in 
the plane of the sky of 22 km s"' (relative to W33 Main), yield- 
ing a total speed of 31 km s"' for W33 B relative to W33 Main. 
In a forthcoming paper, we will compare the gravitational and 
kinetic energies of the clouds to determine if W33 A and W33 B 
are gravitationally bound to W33 Main. 

The massive star forming region G012.88+0.48 appears as- 
sociated with a molecular cloud at a velocity of 30.9 km s"' 
(Fig.fTTI). Thus, both cloud velocity and parallax distance locate 
this source in the Scutum spiral arm. Since previous studies of 



G012.88h-0.48 assumed a distance of 3.6 kpc, its luminosity and 
mass have been overestimated by a factor of two and should be 
revised to L~1.3-10^ L© and M~1.3-10^ M©. 



5. Conclusion 

Trigonometric parallax observations of water masers in the 
G012. 88-1-0. 48 region and in the massive star forming complex 
W33 (containing G012.68-0.18, G012.81-0.19, G012.90-0.24, 
G012.90-0.26) yield distances which are all consistent with 
2.40;^Q j^ kpc. A combination of our distance with CO observa- 
tions of the Galactic plane by iDame et aP (1200 ih allows us to 
locate the W33 complex and G012.88-h0.48 in the Scutum spi- 
ral arm. We show that W33 is a single star forming complex at 
about two-thirds the kinematic distance of 3.7 kpc. Thus, the lu- 
minosity and mass estimates of this region should be lowered to 
L ~ 8-105 L© and M ~ (0.8-8)- 10^ M©. The spectral types in 
the star cluster in W33 Main should be changed by 1.5 points to 
later types, yielding spectral types ranging from 07.5 to Bl .5. 

The luminosity and mass values of G012.88h-0.48 have also 
been overestimated in the past and revised values are L~1.3-105 
L© and M~l. 3-10^ M©. 
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